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ABSTRACT 

In this thorough study we focus on the indirect detection of Dark Matter (DM) 
through the confrontation of unexplained galactic and extragalactic 7 -ray signa¬ 
tures for a low mass DM model. For this, we consider a simple Higgs-portal DM 
model, namely, the inert Higgs doublet model (IHDM) where the Standard Model 
is extended with an additional complex SU(2)i doublet scalar. The stability of 
the DM candidate in this model, i.e., the lightest neutral scalar component of 
the extra doublet, is ensured by imposing discrete Z 2 symmetry. The reduced-y 2 
analysis with the theoretical, experimental and observational constraints suggests 
the best-fit value of DM mass in this model to be ~ 63.5 GeV. We analyse the 
anomalous GeV 7 -ray excess from Galactic Centre in light of the best-fit IHDM 
parameters. We further check the consistency of the best-fit IHDM parameters 
with the Fermi-LAT obtained limits on photon flux for 18 Milky Way dwarf 
spheroidal satellite galaxies (dSphs) known to be mostly dominated by DM. Also 
since the 7 -ray signal from DM annihilation is assumed to be embedded within 
the extragalactic 7 -ray background (EGB), the theoretical calculations of photon 
flux for the best-fit parameter point in the IHDM framework are compared with 
the Fermi-LAT results for diffuse and isotropic EGB for different extragalactic 
and astrophysical background parametrisations. We show that the low mass DM 
in IHDM framework can satisfactorily confront all the observed continuum 7 -ray 
fluxes originated from galactic as well as extragalactic sources. The extensive 
analysis performed in this work is valid for any Higgs-portal model with DM 
mass in the ballpark of that considered in this work. 
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1. Introduction 


The presence of dark matt er (DM) in the universe is now establi s hed following various 


astrophysical and c osmological ( Ade et al. 2014 : Beeeman et ah 1991: iKomatsu et al. 2011 


Massev et ajJ 12007) evidences. The recent data of PLANCK flAde et al. 20141) suggest that 


~ 26.8% of the total mass-energy content of the universe consists of cold (or non-relativistic) 
dark matter whose particle nature is yet to be resolved. The weakly interacting massive 
particle, commonly known as WIMP appears to be the most promising particle candidate 
of cold DM in the universe. 


There are various ongoing experiments for the detection of dark matter both through 
direct and indirect mechanisms. In case of direct detection, the dark matter may scatter 
off a nucleus of a detecting material and in such direct detection experiments attempts 
are made to measure this recoil energy of the nucleus as the signature of dark matter de¬ 
tection. T here are various ongoing direct detection experiments aro und the world such a s 


CDMS II ( Ahmed^jdyiL 2010 : Agnese et ah 2013bl a ), CRESST II (Angloher et al. 


CoGe NT (lAalseth et al. 20111). XENON 100 flAprile et al.1 201 ll. 120121 1. LUX flAkerib et al. 


20121 


20141 1 etc. that use different detection material, e.g. Ge, Si, Xe etc. 


The dark matter in the universe, because of its all pervading nature, may be trapped by 
very massive astrophysical objects such as galactic centre, solar core etc. and may undergo 
multiple scattering with the dense matter present at those sites losing in the process their 
velocity of escape and eventually are trapped inside these bodies. When accumulated in large 
numbers, these trapped dark matter particles may undergo pair annihilation to produce the 
pairs of standard model particles such as qq or ££ as primary or secondary products. Gammas 
can be obtained as secondary products from the pair annihilation of these primary pairs of 
SM fermions (such as hadronisation of bb through 7T°). The indirect dark matter search 
experiments look for these gamma rays or the other SM particles from dark matter pair 
annihilation. 


The satellite borne detectors such as Fermi-LAT or Fermi Large Area Telescope detect 
the gamma rays from the the galactic centre (GC) and inner galaxy regions. Any anomalous 
GeV gamma-ray excess from the Fermi-LAT observation may indicate DM pair-annihilation 
at the galactic centre region in case other known astrophysical phenomena fails to explain 
such excesses. This excess gamma-ray signal can be explained with DM models where DM 
particle annihilates mainly into qq channels with the desired value of thermally averaged 
annihilation cross-section similar to the canonical cross-section of typical thermal production 
of DM. This bump -like feature indicating gamma ray excess is also reported by Fermi- 
LAT collaboration (IMurgial 2014 ) for the gamma rays from the galactic centre region. An 
early analysis of Fermi-LAT data reveals that the gamma rays from the galactic centre 
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region exhibit excesses (bump) in the energy range 


Hooperll2012l;lHooper fc Sla 


wer 


0.3 — 10 GeV flHooper fc Linden! 120 111; 
2013 ; Huang et al. 2013). More involved and modified recent 


analysis flDavlan et a 1.1 12014) including more recent data restricts the range of excess gamma 
to be in the energy region of ~ 1 — 3 GeV. 


The early analyses flHooper fe Linden! 201 ill by Dan Hooper et al. suggest that the 
gamma-ray excesses mentioned above and the morphological features of these excesses can 
be satisfactorily explained by considering DM particles in the mass range o f ~ 30 — 60 GeV 
annihilating only through bb channel. From the latter analysis flHooper fe Slatverl 2013 


Huan g et al. 20131). the authors also gave the best fit value of the dark matter mass to 
be 61 . 8^49 GeV annihilating only into bb pair with thermally averaged annihilation cross- 
section (av) = 3.30lo;4g x 10 26 cm 3 s 1 for explaining the above mentioned gamma ray 
excess. Another analysis flHooperl 201_2) considering the dark matter to have primarily 
annihilated only to rf, yields the dark matter mass in the range ~ 7 — 10 GeV with 
thermally averaged annihila t ion cross-section (c ry) = 5.6 x 10 ^ 27 cm 3 s _1 . A very recent 


analyses flDavlan et al.l 12014; lLaeroix et al. 2014) however indicate that a DM particle with 
mass ~ 31 — 40 GeV and annihilating entirely into bb channel with thermally averaged 
annihilation cross-section (av) = (1.4 — 2.0) x 10 ~ 26 cm 3 s _1 (normalised to a local DM 
density of 0.3 GeV/cm 3 ) can provide much better agreement with the nature of the low 
energy gamma-ray spectrum (with an excess in the energy range ~ 1 — 3 GeV). This new 
analyses also disfavour the possibility of the previous proposition of 10 GeV DM anni¬ 
hilating solely into rr channels. @ There are several attempts (Logan 2011 ; Bucklev et ah 


201 


l3; 


Zhu 


2011 


Marshall fe Primulando 


201 lat Hooper et al.l hold; Bucklev et ah 


20 lit Boucenna fe Profumo 


2013; 


2014; ICahill- Rowlev et al.l2015 


Anchordoaui fe Vlcek 


2011 


Brickiev et al. 


2013 


Dasgupta 


Mo dak et al. 


2013; Guo et al. 


2014; Yu 


Borah fe 


12015 


Ba,nik fc Maiumdar 


2014; 


Okada fc Seto 2014a; 


Ghosh et al. 2015; 


Cheung et al 


Ko et al. 2014; 


2014; 


2014; 


2014; 


Basak fc Mondal 


3alzs fe Li 12014; Agrawal et ah 2014b,c; Izaguirre et ah 


.2014a 


2014; iBerlin et al~ 


2014a: 


2014 


2014 


Cerdeo et a' 


Arina et al 


2015; 


Ipek et al. 


Huang et al. 


b; Alves et 


2014 


. 3oehm et a. 


2014; 


Wang fc Han 


2014 


Ko fc Tang 20151: Agrawal et al.l2014a; 


Fields et al. 


Okada. fc Seto 


2014b 


Boehm ct al 


20141 


al 


Berlin ct al. 


2014b 


Abdullah et al. 


2014; 


Martin ct al. 


2014 


Cline et al 


2014: 


Detmold et al 


2014- 


Chang fc Ng 


2014; 


Bell et ah 


2015 


Cao et ah 


iHeikinheimo fc Spethmann 2014: Agashe et al. 12014: Ghorbani fc Ghorbani 


2014 ; lFrevtsis et al 


2015 


2014 ; ICerdeno et al. 20151) to propo se DM models and studying various a spects. In a more 


recent analyses icalore et al. 2014a b; Murgia 2014 ; Agrawal et al. 2014a ) the galactic cen¬ 
tre excess has been reanalysed considering several distinct galactic diffuse emission models 


Cosmic ray positron data also disfavour the possibility of DM annihilating into r lepton pairs with the 
proposed mass of DM (~ 10 GeV) and annihilation cross-section to accommodate GeV gamma-ray excess. 
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and the allowed DM mass range for the generation of such galactic centre gamma ray excess 
is severely relaxed. The pr eferred mass range of DM anni hilating solely to bb channel is 
derived to be 35 — 165 GeV (Agrawal et ah 2014a). In Ref. f Calore et al. 2014b) the Fermi- 
LAT galactic centre GeV excess is interpreted with DM mass allowed up to 74 GeV with 
bb annihilation channel. Alternative propositions other than annihilating DM explanations 
such as unresolved millisecond pulsars to be the main origin of this observed gamma-ray 
excess have been discarded since the observed anomalous gamma-ray emission extends much 
beyond the central stellar cluster. 

Besides galactic centre, the dwarf spheroidal galaxies or dSphs also are very rich in dark 
matter. These are faint companion galaxies of Milky Way. The mass to luminosity ratio (^) 


are found to be much higher than 


M 

L 


0 


where 


M 

L 


0 


denotes the mass to luminosity ratio of 


the sun indicating that these are rich in dark matter. These dark matter can pain annihilate 


analyses ( 

Abdo et al. 

2010 a: 

Ackermann et al. 

2011 ; 

1 .7 - ^ 7 — -- • 

Gerineer-Sametli & Koushiappas 

2011 ; 

Mazziotta et al. 

2012 

Gerineer-Sameth & KoushiapDas 2012; 

Ackermann et al. 

2014 

) per- 


formed on several dSphs to constrain DM annihilation. 

Apart from the galactic cases, the observed y-ray signal by Fermi-LAT from the ex- 
tragalactic source s also may conta i n the signature of the dark matter ann i hilation at ex¬ 
tragalactic sites (Ullio et al. 2002 ; Bergstrom et al. 2001 : Gao et al. 1991 : Stecker 1978 


Taylor & Silk 2003 


Ng et al. 


2 0141 ). The signal may also have embedded in it the y- 


tempts ( 

Brinemann et al. 

2014; 

Cholis et al. 

2014; 

Tavakoli e 

al. 

2014; 

1 

Sefusat 

ti et al 

2014; 

Aiello et al. 

2015; 

Di Mauro & Donato 

2015 

Di Maurol 

2015; 

Ackermann et al. 

2015a) 

to ex- 


tract DM signal from such extragalactic 7 -ray background (EGB) and to provide limits on 
DM annihilation cross sections for different DM masses. This requires proper modelling 
of the extragalactic parameters as well as proper knowledge about the other astrophysical 
backgrounds that contribute dominantly to the EGB signal. With the analyses of new data 
collected by Fermi-LAT mission, a much detailed and clear picture of EGB has been put 
forward. The information regarding the astrophysical sources such as BL Lacs, millisecond 
pulsars, star forming galaxies, radio galaxies etc. which possibly contribute to this EGB are 
unveiled from various observations in radio and gamma wavelengths. As Fermi-LAT collects 
more data, one can precisely measure the EGB spectra and put stringent constraints on DM 
annihilation cross section. This constraints are contemporary to that obtained from dSphs 
and galactic centre regions and may, in principle, put limits on various DM models in future. 


A number of particle physics models for the dark matter candidate has been proposed 
and studied in literature. They include various extensions of Standard Model flCheng et al. 
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2002: Servant & Ta.it 

2003; 

Duffv & van Bibber 

2009; 

Randall & Sundrum 1999: Ma 

2006 

) 

whose DM phenomenologies 

(Modak & Maiunifh 

11 201 

3k Modak 20141 are explored at length. 

Amongst them, the Higgs-portal models such as singlet scalar DM (Silveira & Zee 

1985 

; 

5 

) 

Veltman & Yndurain 

1989; 

Burgess et al. 2001; I 

Barger et al 

20081: Gonderinger et al. 

2010 ) 

inert Higgs doublet model HHDM) (Deshpande & Ma 

19781 

Lopez Honorez & Yaguna 

2011 


and two Higgs doublet model (2HDM) (IBranco et al 


20121) , singlet vec tor DM flKanemura et al. 


2010 ; Lebedev et al. 20121; Diouadi et al. 2012 b singlet fermionic DM ( Kim et al. 2008 1 could 
be of particular interest for the present scenario in explaining the observed anomalous gamma 
emission by Fermi-LAT. The Higgs-portal models are interesting to study since the low mass 
DM candidates of these models annihilate into quark pairs with the cross-section in the right 
ball park of thermal production. A special feature of these types of models is that the DM 
candidates in these models exhibit resonance phenomena when their masses reach the value 
of ~ half of the Higgs mass while satisfying bounds given by PLANCK experiment (relic 
density) and dark matter direct detection experiments. 

In this wor k, we focus on t he I nert Higgs Doublet model (IHDM), proposed by Desh- 


pande and Ma flDeshpande fe Malll978f ) and confront the recently observed gamma-ray ex¬ 


cesses from galactic centre region with the dark matter candidate in this model. We also 
explore the possibilities of the observation of gamma rays from 25 dwarf spheroidal galax¬ 
ies by this IHDM dark matter candidate. I11 addition, we study the extragalactic gamma 
ray signals obtained by Fermi-LAT with this inert Higgs doublet dark matter. In the in¬ 
ert Higgs doublet model or IHDM an extra scalar doublet is added to the Standard Model 
(SM) which is assumed to develope zero vacuum expectation value after spontaneous sym- 


Krawczvk et al. bold: 


noil 


DM p 

lenomenology 

(G01 

rdelis et al. 

2013 

: Gustafsson et al. 

2012 

: L 

ooez Honorez & Yaguna 

201 ll: 

Ha.mbve & Tvtgat 

2008 

; P 

Lgrawal et al. 

2009; 

Harnbve et al. 

2009: Nezri et al. 

2009; 

Andreas et al. 

2009; 

Arina et al 

20091 I 

^ooes 

Honorez & Yagunal 

2010 ; 

Melfo et al. 

2011 ; 


Gustafsson et a 


20121; Swiezewska. fe Krawczvkl 2013bl . 


a]; 


20121 


Dolle fe Su 


Lundstrom et al 


Kanemura et al 


20091: 


Cao et al. 


2011 


200 ?]; 


Arhrib et al. 


Dolle et al. 


2010|). For the dark matter candidate in this IHDM framework, a reduced y 2 analysis is 
performed considering all the above data and constraints and the best fit values for dark 
matter mass, annihilation cross-section and other parameters of the model (various coupling 
constants) are found out. In the present work we adopt those best fit values as the bench¬ 
mark point and study the Fermi-LAT gamma-ray flux results from both galactic (galactic 
centre, dSphs) and extragalactic sources. 


The paper is organised as follows. In Section [2] the theoretical framework of the IHDM is 
briefly described. Also the theoretical, observational and experimental constraints imposed 
on this model are discussed in this section. Confronting the observed gamma-ray excess from 
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galactic centre in this model framework with a detailed study of the computed gamma-ray 
spectra is performed in Section [3j We compare the calculated results with the bin-by-bin 
upper limits on photon energy spectra for various Milky Way dSphs in Section SI The Sec¬ 
tion ED contains the confrontation of the extragalactic 7 -ray background with calculated pho¬ 
ton spectra in IHDM considering different extragalactic parametrisations and astrophysical 
non-DM backgrounds. In Section [ 6 ] we summarise our study and some important conclusions 
have been drawn. 


2. Inert Higgs Doublet Model (IHDM) framework 


The Inert Higgs Doublet Model is one of the simplest extensions of Standard Model 
(SM) Higgs sector where an additional complex scalar doublet, <f>, odd under the discrete 
symmetry Z 2 , is considered alongwith the SM Higgs doublet, H\. After spontaneous sym¬ 
metry breaking, while the SM Higgs gets a vacuum expectation value (vev) v, the additional 
doublet does not acquire any vev. Thus under Z 2 symmetry $ —» — $ and H —>■ H (even 
under Z 2 ) and after symmetry breaking the two doublets H and $ can be expanded as, 


H = 


1 

V2 


G + 

(v + h° + iG°) 


$ = 


\/2 


H+ 

(H° + L4°) 


( 2 . 1 ) 


where Cr and G° are charged and neutral Goldstone bosons respectively. Note that vev of 
these scalar doublet fields are (H) = v/y/2 (v — 246 GeV) and ($) = 0. 

With the unbroken Z 2 symmetry the model has a CP-even neutral scalar H°, a CP-odd 
neutral scalar A 0 , and a pair of charged scalars H ± . Since the Z 2 symmetry excludes the 
couplings of fermions with H°, A 0 , H ± , the decay of the latter particles to fermions are thus 
prevented. This ensures the stability of lightest neutral scalar (. H° or A 0 ) and hence the 
lightest among these two can serve as a possible DM candidate. Either H° or A 0 is chosen as 
the lightest inert particle or LIP and is the candidate of dark matter in the present model. 

The most general tree-level scalar potential of IHDM consistent with imposed Z 2 sym¬ 
metry can be written as, 


Ho = /i?|i/| 2 + / r2|$|2 + A 1 |Lr| 4 + A 2 |<I)| 4 + A3|i/| 2 |4>| 2 + A4|i/ t <h| 2 + ^ 


(i / t $) 2 + h.c. , ( 2 . 2 ) 


where /i*s and AjS denote various coupling parameters. The model has set of six parameters, 
namely 


{M h 0, M h 0, M a 0, M h ±, X l , A 2 }, 


( 2 . 3 ) 
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where M h o, M H o, M^o, M H ± are the masses of Higgs h, CP-even scalar H°, pseudo-scalar 
A 0 and charged scalars H ± . X L and A 5 are in couplings given by, 


g (^ 3 + -^4 + A 5 ), 
As = - (^ 3 + -^4 — A5 ). 


(2.4) 

(2.5) 


The parameters A^ or As denote the coupling strength for H°H°h° (if H° is considered to 
be the lightest inert particle or LIP) or A°A°h° (if A 0 is the LIP). 


A detailed study of this model had been done in Ref. fiArhrib et al.l 20141) where the 
authors made use of the various constraints available from DM experiments and other re¬ 
sults and made a y 2 analysis with the IHDM theory of the dark matter mentioned above. 
Considering all possible experimental and theoretical constraints such as Planck limits, di¬ 
rect detection constraints, unitarity, perturbativity, etc. as also constraints from LHC, they 
provide the best fit values of the quantities M h 0 , M H ±, M H 0 , M A 0 and the parameters A l, 
X 2 (shown in Table [T]). 


3. Confronting Gamma Ray flux from Galactic Centre in this framework 


The differential gamma-ray flux due to the annihilating DM coming from the galactic 
DM halo per unit solid angle can be written as flCirelli et al. [20 11|), 


d<h 7 1 (crv) f dNf 2 

d{ldE 1 M| O A o d £ 7 r ° Po ’ 


where 0 is the mass of DM candidate H° and a = 1 for the present DM candidate (self- 
conjugated). In Eq. 13.11 r 0 and p Q are the distance of the solar system from the galactic 
centre and the local DM halo density respectively. The factor J in Eq. 13.11 gives the total 
dark matter content at the target and is given by, 


J 


ds_ /p(r)A 2 
r& V Pe J 


(3.2) 


where p Q and p(r) are respectively the DM density at solar region and the density at a radial 
distance r from GC and r is expressed in terms of line of sight, s as, 


(s 2 + Tq — 2 sr 0 cos l cos 6 ) ^ 
(s 2 + Tq — 2 sr Q cos 6) 1 ^ 


(galactic l, b coordinate) 
(galactic r, 6 coordinate) 


(3.3) 














In the above p(r) is the DM h alo pro file and f or a g eneralised NFW DM halo the analytical 
expression can be given (INavarro et al.lIl99 6l 119971) as, 


p(r) = 


Pe 


[r/r c p[l + (r/r c ) 7 ]l 3 - 0-7 ) 


(3.4) 


and Pq ~ 0.4 GeV/cm 3 . In the above relation, the values of the parameters r c a nd y are taken 


to be 20 kpc and 1.26 respe ctive ly in the p resent calculation following Ref. flDavlan et ah 


20141) . We use micrOMEGAs ( Belanger et al.l 1201 ll, 120141 ) code to compute various DM ob¬ 


servables in this model. 


The calculated branching ratio for the channel LIP LIP —» bb is found to be 69. 

The branching ratios for other annihilation channels in case of the present LIP dark matter 
are also computed and they are tabulated in Table [21 Since the DM mass is close to half 
of the SM Higgs (~ rrih/2) like particle discovered by LHC, we will have resonance effect 
in obtaining the required cross-section (2.37 x 10 -26 cm 3 s _1 ) with bb as dominant channel 
for the typical thermal production of DM. Also DM-nucleon scattering cross-section for 
the chosen benchmark point in this model is about 8.89 x 10” 11 pb (averaged value per 
nucleon for interaction with Xe nucleus) which is just under the the present bounds for 
XENON 100 and LUX exper iments. This is shown in Fig. [H The future DM direct searches 
like XENON1T flAprilel 120131 ) can easily probe this point as seen from Fig. Q] 


The Fermi-LAT data for gam ma-ray flux from the inner 5° surrounding of the galactic 
centre have been studied in Ref. (lHooperll2012n . The known y-ray sources in this region can 
be fo und in (extracted from) Fermi Second Source Catalog (2FGL) ([Fermi-LAT Collaboration 
2 012l) . Although Fermi Third Source Catalog (3FGL) has recently been made available 
(Fermi- LA T Collaboration 120151) . but no analysis of the backgroimd or analysis for y-ray 
from other known processes at the region of interest has been reported yet. Also cosmic ray 
interaction with gas distributed in this galactic region produces neutral pions that subse¬ 
quently decay to produce enormous amount of y-ray. This is a viable mechanism for known 
disc template emission. Now, the spectral and morphological feature of the photon flux from 
inner 5° subtending the GC after subtracting the contribution from both the known sources 
of Fermi Second Source Catalog and disc template emission shows a ‘bumpy’ nature in y-ray 
energy ranging from ~ 300 MeV to ~ 10 GeV. The count drops significantly after 10 GeV 
of y-ray energy. 


We have computed the y-ray spectrum from inner 5° subtending the GC from DM 
annihilation within the present framework of inert Higgs doublet model for dark matter 
particle. The chosen benchmark points for the parameters of the model such as dark matter 
mass, coupling constants etc. are given in Table [Q The flux have been computed for the 
generalised NFW DM halo profile. Also included in the calculation are the contributions from 
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both point source and galactic ridge emission. The total calculated flux is then compared 
with the observed residual photon flux and the results are shown in Fig. [2] In Fig. [2] the 
green-coloured and blue-coloured lines represent the fluxes for point source and galactic ridge 
emission respectively. The y-ray spectrum for DM annihilation for the benchmark points 
mentioned above, are shown in purple line whereas the black line is for the total y-ray flux 
obtained by summing over all the fluxes represented by green, blue and purple coloured lines 
in Fig. [2] For these calculations we have adopted the generalised NFW (gNFW) halo profile 
with 7 = 1.26. The total flux (black line) is then compared with the observed residual 
emission data. These observed data are denoted by red-coloured points in this figure. It is 
clear from Fig. [2] that our computation of total residual gamma emission (black line) agrees 
satisfactorily with the observational results. 


In a recent analysis of y-ray flux from GC region where the analysis is only for the GC 
gamma r ay (subtrac ting all possible contributions from other known astrophysical sources), 
flDavlan et al 120141) an excess of gamma ray in the gamma energy region of ~ 1 — 3 GeV 
has been reported. This excess is shown in the left panel of Fig. [3] with the red-coloured 
points. It is suggested in the same analysis that in order to explain this anomalous gamma 
ray excess from dark matter annihilation scenario, the DM mass should be in the range of 
31-40 GeV which is to annihilate purely into bb pair. We calculated these gamma ray fluxes 
in our framework of inert Higgs doublet dark matter for the dark matter mass of ~ 63.5 
GeV (adopted from Table [1] (benchmark point)) and compared our results with the analysed 
data points mentioned above (red coloured points shown in the left panel of Fig. [3]). In 
the left panel of Fig. |31 the green coloured line represents the present calculation. These 
calculations are performed considering gNFW halo profile with the halo parameter y = 1.26. 
It is to be noted from the left panel of Fig. [3] that although the morphological feature of 
the spectrum from our calculation (green line) is similar in nature to that obtained from 
the analysis of Fermi-LAT observational data (red-coloured points) from GC, the position 
of the maxima of excess gamma ray in our calculation is shifted to somewhat higher energy 
at ~ 3.1 GeV instead of being within the expected energy range of ~ 1 — 3 GeV as obtained 
from the Fermi-LAT data. However the calculated position of the peak (green line) is at 
~ 2.84 x 1CV 6 GeV/cm 2 /s/sr which is in the same ballpark of the observed peak (red points). 


More detailed analysis of the observed y-ray excess reveals (JD a y l an et ah 2014!) a further 
anomaly between the gamma ray spectra for the gamma rays from galactic East-West region 
and from galactic North-South region. The North-South region is designated as \b\ < |h|, 
where b and £ are the galactic latitude and longitude respectively and for the East-West 
region, |6| > \£\. The two spectra are shown in the right panel of Fig. El The red coloured 
points are for gamma spectrum from galactic East-West region while the blue coloured points 
represent the spectrum from North-South region. As can be seen from the right panel of 
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Fig. [3] the gamma flux from the present calculation (shown by green line in the right panel 
of Fig. |3]) agrees more satisfactorily with the “North-South” gamma emission spectrum than 
the “East-West” spectrum. 

The systematic uncertainty for the estimation of background model pro vided by the 
Fermi-LAT is very large compared to the statistical uncertainty. Attempts flCalore et al. 
2014bl ; iMurgia 2014) are made to quantify such systematics for galactic central region. We 
confront the flux obtained for our IHDM benchmark scenario (Table ED) using these two 
complementary approaches. We adopt the proposed method of Ref. flAgrawal et al.ll2014ai) 
for the uncertainty estimation of J-factor. This involves estimation of halo profile uncertainty 
and the uncertainty in local DM density. The density profile index 7 , is estimated to be 
7 = 1.2 ± 0.1 from different galactic diffuse emission models and p Q is estimated to be 
p 0 = 0.4 ± 0.2 GeV/cm 3 for different normalisation of halo.Defining the J-factor as 


J = 


AD 


'AO 


J('0)c?D — I X Jcanonical ) 


(3.5) 


where ./canonical is the central value of J and the factor JJ signifies the deviation from the 
canonical halo profile due to the uncertainties of the profile. In the above AD is the region 
of interest (ROI) for a given analysis. The astrophysical factor J(V’) is same as that in 
Ea.IiOl 


I 11 Ref. flCalore et ah 20 14b| ) a thoroughly analysis of Fermi-LAT data in the inner galaxy 
region over the photon energy ranging from 300 MeV - 500 GeV has been made where the 
chosen region of interest (ROI) is extended to a 20° x 20° (—20° < l < 20°, —20° < b < 20°) 
square region surrounding the galactic centre. The inner galactic latitude of 2° (| 6 | < 2°) 
has been masked out. For the canonical halo profile with r 0 = 8.5 kpc, p Q = 0.4 GeV/cm 3 , 
r c = 20 kpc and 7 = 1 . 20 . The numerical value of J ca nonicai is found out to be 2.0 x 10 23 
GeV 2 cm -5 for the canonical halo. The uncertainty (J G [0.19,3]) in the J factor is found 
out to be in their analysis. 

Plugging in all of the canonical parameters for the dark matter halo model, we compute 
the canonical 7 -ray flux obtained from annihilation channels of ~ 63.5 GeV LIP DM (adopted 
from Table [2] (benchmark point)). The resulting plot is shown in the left panel of Fig. [4] 
by black coloured line. The black line is obtained considering only the canonical value 
J G [0.19,3] of the J-factor. The red coloured points denote the residual spectrum of 7 -ray 


excess in the galactic centre with highly correlated errors obtained from Ref. flCalore et al. 


2014bh . We repeat the calculations by also taking the uncertainties J G [0.19,3] in the J 
factor. The blue and the green coloured lines in the left panel of Fig. [4] denote the galactic 
7 -ray flux from ~ 63.5 GeV DM annihilation for the maximum value of the uncertainty 
Jmax ~ 3.0 and the minimum value of the uncertainty J min ~ 0.19 respectively. From the 
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left panel of Fig. [4] one sees that the uncertainty factor J needs to be smaller than unity to 
have a better fit to the data. 


On the other hand, recently the Fermi collaboration (Murgia 2014 ) has also studied the 
region surrounding the galactic centre using different background models of galactic diffuse 
emission. The fit to the galactic centre y-ray data is observed to improve very significantly 
when an additional contribution similar to that from dark matter annihilation is added. The 
Fermi collaboration has chosen a di fferent reg ion (15° x 15°) surrounding the galactic centre 
smaller than that chosen by Ref. (ICalore et al.ll2014b ). However for their analysis, unlike in 
Ref. (ICalore et al. 2014b), in this case, the galactic centre is not masked out. Based on their 
preliminary analysis the Fermi collaboration has reported four best fit y-ray spectra for the 
four distinct choices of background models for galactic diffuse emission. The nature of these 
four best fit y-ray spectra differ notably after a few GeV photon energy. The obtained y-ray 
spectra is found to yield much more conservative measurement of the systematic uncertainty. 
Although the Fermi h as ana l ysed the data using NFW halo profiles with slope values of 1.0 
and 1.2, Ref. (lAgrawal et al.l I2014aj ) has chosen a more conservative approach on y factor 
and set it to 1.2 ± 0.1. Also the value of p Q is chosen to be 0.4 ± 0.2 GeV/cm 3 . Keeping the 
parameters of the J-factor fixed, one would obtain in this case, J ca nomcai = 1-58 x 10 24 GeV 2 
The uncertainty in the J factor as obtained from such parametrisation is J G [0.14,4.0] 


cm 


-5 


for the Fermi analysis. 


We make an estimation of the galactic y-ray excess from the annihilation of dark matter 
as chosen in this work (benchmark point in Table [T]) and confront the results with the y- 
ray spectra obtained by Fermi collaboration. We make the comparison of the calculated 
spectrum with that reported by Fermi after preliminary analysis and show it in the right 
panel of Fig. [4j The black line in the right panel of Fig. [4] denotes the photon flux obtained 
using the canonical parametrisation of the halo profile (with J = 1). The blue and the green 
coloured lines in the right panel of Fig. [I] represent the galactic y-ray flux calculated for the 
maximum uncertainty J max ~ 4.0 and the minimum uncertainty J m ; n ~ 0.14 respectively. 
Also shown in the right panel of Fig. [4] the Fermi analysis results with its upper and lower 
bounds. § From the right panel of Fig. [4] it appears that unlike the previous analysis (with 
—20° < £,b < 20°) the uncertainty factor J should be more than unity in order to provide a 
better fit to the data for the considered IHDM benchmark point. 


§Out of the four spectra given by the preliminary analysis of Fermi, we choose only a particular spectrum 
which provides the best fittings for dark matters with low masses since our interest in this paper is on the 
low mass region in IHDM. 
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4. Confronting Gamma-ray flux from Dwarf Spheroidal Galaxies in this 

framework 


One of the most promising targets for the search of dark matter via indirect detection 
( 7 -ray) is the dwarf spheroidal galaxies (dSphs) of the Milky Way. The dSphs are consid¬ 
ered to be promising for the study of DM phenomenology because of their proximity, low 
astrophysical background and huge amount of DM content. 


The satellite-bourne gamma-ray experiment, Fermi-LAT, search for the y-ray sky in the 

energy range spanning from ~ 

500 MeV - 500 GeV ( 

Atwood et al. 

2009 b In a more recent 

study by Ackermann et al. ( 

Ackermann et al. 

20 D 

1), 4-year gamma ray data of Fermi- 


LAT on dSphs (04-08-2008 to 04-08-2012) with energy ranging from 500 MeV to 500 GeV 
have been chosen for studying 25 independent Milky Way dSphs galaxies. The chosen dSphs 
galaxies are Bootes I, Bootes II, Bootes III, Canes Venatici I, Canes Venatici II, Canis Major, 
Carina, Coma Berenices, Draco, Fornax, Hercules, Leo I, Leo II, Leo IV, Leo V, Pisces II, 
Sagittarius, Sculptor, Segue 1, Segue 2, Sextans, Ursa Major I, Ursa Major II, Ursa Minor 
and Willman 1. The galactic coordinates as well as the radial distances from the galactic 
centre of these dwarf galaxies are tabulated in Table [3j From the analysis of their data, they 
set robust upper limits on DM annihilation cross-section for different DM masses. In giving 
this limits, they consider the DM pair annihilate predominantly to bb as also rr and other 
channels. 


It is possible to assess the total DM content of dSphs galaxies from the dynami¬ 
cal modeling of the stellar density of the dwar f galaxies and th e velocity dispersion pro- 
flies ( Battaglia et al. 2013 ; Walker et ah 2009b ; Wolf et ah 2010 b The dynamical masses 
of these dwarf galaxies are measured only from stellar velocity dispersion and half-light ra¬ 
dius. The calculated total mass within the half-light radius for a dSphs galaxy is used to 
obtain the integrated J-factor of that dSphs galaxy. Both the total mass within the half- 
light rad ius and the ./-factor are found to be almost indepen dent on the choice of DM halo 
profiles ( Martinez et ah 2009 ; Strigari 2013 ; Martinez 2013 b Out of the 25 independent 
dSphs me ntio ne d ear lier, J-factors of only 18 dSphs are determined using stellar kinematics 
data (IMart inezI 120131 ) while other seven lack proper statistical significances. Using these 
uncertainties of J-factors, the upper bounds on DM annihilation cross-section for various 
DM masses have been derived with 95% CL. 


We calculate the 7 -ray flux for all the dwarf galaxies. The 7 -ray spectrum, ^ can be 
obtained for a given DM mass. The different particle processes for the calculation of ^ are 
tabulated in Table [2j We compute ^ for our benchmark scenario and using the integrated 
J-factor for a particular dSph, the maximum value of the velocity-averaged annihilation 
cross-section ((cru) max ) is estimated from the upper bound of the flux (LHS of Eq. 13.IK of 
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that dSph. In this way the upper bounds of annihilation cross-section are computed for all 
the 18 dSphs considered and they are tabulated in Table [3j 


We now compute the 7 -ray flux for all the 18 dSphs considered using Eq. 13.11 The 
maximum, minimum and central values of integrated J-factor for each dwarf g alaxy which 
are measured from stellar kinematics data are tabulated in Table l3l (jAckermann et al. 2014 b 


The integration to find J-factor involves integration over a solid angle All of ~ 2.4 x 10 " 4 sr 
(the field of view of Fermi-LAT is within the angular radius of 0.5°). The results are shown 
in the plots given in Fig. E3We also show in Fig. ED the experimentally obtained bin-by-bin 
upper limits of the gamma ray energy flux at 95% CL from each dwarf galaxy by downward 
red-coloured arrows. The flux of a dSph are compared with the upper bound of the flux in 
each energy bin. The spread (band) of this flux shown by green indicate the upper and lower 
limits of the flux when calculated with the upper and lower limits of integrated J-factors. 
The photon flux calculated using the central value of J-factors are shown by blue lines in 
these figures. 


5. Confronting Extragalactic gamma ray background in this framework 

The 7 -rays can also be emitted from DM annihilation i n extragalactic sources and such t 


rays 


can be probed for extragalactic DM and their origins flUllio et al 


2001 


2005 


Gao et al. 1991 


2002I; 

3erestrom et al. 

; Ando 

2005: Oda et al. 


Pi eri et al . I 2 OO 8 I. Such gamma rays from extragalactic sources of DM can remain hid¬ 


den in the huge background of the o bservationally meas ured gamm a flux by satellite-born e 
experiments such as SAS-2 satellite (Fichtel et al. 1978 1. EGRET f Sreekumar et ah 1998 1. 
Fermi-LAT flAbdo et al.l l2010bl: I Ackerman n et aLll2015bfl . In order to extract information 
regarding the extragalactic signature of gamma rays from the background one should be 
able to understand and subtract the galactic astrophysical components, other sources that 
may contribute to the background and the backgrounds that the detector may give rise to, 
in the process of detection. After this process of subtraction the residual gamma-ray signal 
thus obtained is found to be diffuse and isotropic in nature and is known as diffuse isotropic 
gamma-ray background (DIGRB). Recently in the light of 50-month Fermi-LAT data an up¬ 
dated tight constraint on DM annihilation is given with the modelling of integrated emission 
of blazars with such diffuse background absorption flAiello et al] 120151 1. This may also be 
noted that the DIGRB thus obtained embeds in it the irreducible contributions from galactic 
origin as well. In this section we estimate such diffuse isotropic gamma ray background or 
DIGRB for the case of dark matter annihilating into gamma rays in the framework of the 
chosen IHDM dark matter candidate. We then compare our theoretical calculations with 
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EGRET and Fermi-LAT results for extragalactic DIGRB. 


5.1. Formalism 


The number of photons which are isotropically emitted from the volume element dV 
in time interval dt, in energy range dE and are collected by detector with effective area dA 
during time interval dto with redshifted energy range dEo can be given by flUllio et al.l 2002), 


dN 7 = e~ T ^ Eo) 


(1 + zf 


dM 4^(M, z) ^ (E, M, z) 


dM 


dE 


dV dA 


An[R 0 S k (r)}‘- 


dE 0 dt 0 . (5.1) 


In the above, the volume element dV at a redshift z is given as 

[R„S h (r)] 2 R 0 , 

d\ — 7— T7 - dr ml detector , 

(1 + zy 


(5.2) 


where dD dete ctor = sin ddddcj) denotes the angular acceptance of the detector. In Eq. 15.11 
and Eq. 15.21 S k (r) is given by the Robertson Walker metric for homogeneous and isotropic 


universe 


ds 2 = c 2 dt 2 — R 2 (t) [dr 2 + Sl(r)(dd 2 + sin 2 0d(j) 2 )] , 
where S k (r) signifies the spatial curvature of the universe. In Eq. 15.11 d ' N ~ 


(5.3) 


X 

dE 


(E, M, z) is the 

differential photon energy spectrum for a generic halo with mass M at some redshift z. The 
term ^g(M, z) is the halo mass function and is defined as number density of bound objects 
with mass M at redshift z. The term e -T ( z ’ £ '°) represents the attenuation of extragalactic 
y-rays which may come from the absorption of these high energy y-rays on the extragalactic 
background light (EBL) and t(z, Eq) is the optical depth which is function of z and Eq. 
The energy and redshift dependence of this attenuation factor is shown in Fig. [6] Detailed 


studies regarding this attenuation are given in Ref. flCirelli et a 


backg round we have adopted the minimal model (Dominguez et al 


2011 


2011^ For ultraviolet 


Franceschini et. al. 


20081 ) obtained after a recent study on blazars. Note that dt 0 dE 0 in Eq. 15.11 is given by 
dtdE = ^ z y (l+z)dE 0 , where t 0 and E 0 are the time and energy respectively corresponding 
to the redshift z = 0. Summing over all the above contributions, the diffuse extragalactic 
y-ray flux due to DM annihilation, can be written as, 


dcj). 


7 


dN~ 


dEo dA dD dto dEo 


= — J dr R 0 e~ T{z ’ Eo) 

c r e ~r{z,E 0 ) f 

An J ^ Hq h{z) 


dM -$t(M, z) (E 0 (1 + z), M, z) 


dM 


dE 


dM ^(M, z) ^ (Eq (1 + z), M, z) ((3.4) 


dM 


dE 
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where c is the speed of light in vacuum, Hq is the Hubble constant at the present epoch and 
h(z) is written as (for spatially flat universe (fl*. = 0 )) 


h(z ) — -\/njvr(l + z) 3 + Ha- 


(5.5) 


where Qm and 11 a are respectively the matter and dark energy densities normalised to the 
critical density of the universe. Th e cos mol ogical da rk matter halo function ^( M,z ) in 
Eq. 15.41 can be written in the form flPress fe Schechterlll974l ). 


dn 


Po,r 


dM M 2 




d log v 
dlogM 


(5.6) 


where po,m is the comoving background matter density. In the above, the parameter v = 
8 sc /a(M), defined as the ratio between the critical overdensity for spherical collapse <5 SC 
(~ 1.686) and cr(M) denotes the variance or the root mean square density fluctuations of 
the linear density held in sphere that contains the mean mass M. The term a 2 (M) can be 
writt en in terms of the linear power spectrum P{k) of the fluctuations flSheth fe Tormen 
1999 1 as, 

cr 2 (M) = J d 3 kW 2 {kR)P(k) , (5.7) 

where W ( kR ) is the Fourier transform of the top hat window function and R is the comoving 
length scale. For collapsed halos, the mass is found to be in the form M ~ ( 4 / 3 ) 7 r R 3 p c (z c ) 
with z c being the redshift where collapsing of halos occurs. The power spectrum P(k ) can be 
parametrised as P{k) oc k n T 2 (k ) where n is the spectral index and T is the transfer function 
that incorporates the effect of scale dependency of the primordial power spectrum generated 
during inflation. This transfer function depends on the nature of DM and baryon density 
in the universe. Thus the transfer function can be calculated from the cosmic microwave 
background data. The variation of the power spectrum P{k ) with wavenumber k for different 
redshifts is shown in the left panel of Fig. 0 The right panel of Fig. 0 corresponds to the 
plot showing the variation of variance a with halo mass M for different values of redshift. 
The function /(z/) in Eq. 15.61 k nown as the multiplicity function, can be modelled in the 
ellipsoidal collapse model flSheth & T orm en 1 9991 ) by, 


vfty) =2A[1 + 


u /2p 


~2n ) 


1/2 


exp — 


v 


,/2 


(5.8) 


where v' = y/au, a = 0.707, p = 0.3 are o btained by fitting Eq. 15.61 to IV-body simulation 
of Virgo consortium (Jenkins et ah Il998l h The value of A is obtained to be 0.3222 For 
the choice of parameter val ues, a = l. p = 0 and A = 0.5, Eq. 15.61 reduces to the original 
Press-Schechter theory flPress fe Schechterl |l974j). It is found in IV-body simulations that 
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the estimations of higher and lower mass halos differ from that predicted by Press-Schechter 
model. This problem can be handled in Sheth-Torman model by considering ellipsoidal 
collapse model instead of spherical one. 

In the left panel of Fig. [ 8 ] the variations of the fraction of mass collapsed or /(er) in 
the ellipsoidal collapse model with redshifts z and the halo mass M are shown. Note that 
/(er) as shown in the left panel of Fig. [ 8 ] can be obtained by simple transformation of f(u) 


2 \P i 


— exp 


' 2<t 2 


In 


by plugging in u = S sc /a(M ) and is given by /(er) = A^J ^ 
the right panel of Fig. [ 8 ] we have s hown t h e variations of the considered halo mass function 
dn/dM of Sheth-Torman model flSheth & To r rne n 119991 ) with redshift z as well as with the 
halo mas s M. All the numer ical calculations related to Fig. [ 8 ] have been performed using 
HMFcalc ( Murray et al. 2013) code. 


For the halo profile we have chosen NFW halo profile (INavarro et ah 19961 19971 ) given 


by, 


P(r) = p s g(r/r s ) = p s 


2 ’ 


(5.9) 


(5.10) 


^f + f 

Any DM halo of mass Mh enclosed at a radius 77 is, 

M h = 4 tt p s rlf{r s /r h ) , 
where f(x) = x 3 [ln(l + x~ l ) — (1 + x) -1 ]. 

Also a halo of mass M at some redshift z can be written in terms of mean background 
p(z) as, 

A'TT 

rp{z)R vir , (5.H) 


,, 4tt . 


where R mr is the virial radius defined as the radius within which the total halo mass M 
is contained with mean halo density A m T p{z). The term A V i r is the virial overdensity with 
respect to the mean matter density which may depend on the cosmological parameters but 
indepen dent of halo ma ss M/, . For the flat cosmology, A„j r (z) can be cast into the following 
form (1 Bryan fe Norman 19981 ). 


A vir ~ (187T 2 + 82 d — 39 d 2 ) , 

with d = d(z) = 7 n % 4 vvrT 7 T- t — 1. We choose the value of A, 


(5.12) 


(Q m (l+2:) 3 +r2A) 


The 7 -ray energy spectrum c -^- (E 0 (1 + z),M,z) (Eq. 15.4[) for the gamma-ray emitted 
inside a halo of mass M at redshift z is written to the form, 


dAfy 

He 


(E, M, z) — 


(av) dNy(E) 


dCyi r V{c[ nr ) [jf] I d 3 r g 2 (r/a ) , 


2 dE 


(5.13) 
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where (av) is the thermally averaged value of annihilation cross-section times the relative ve¬ 
locity, dA ( ]l R - is the differential 7 -ray energy spectrum produced per unit annihilation of dark 
matter and M x is the mass of dark matter. The log-normal distribution V(c , nir ) of th e concen ¬ 
trati on parameter c V i r around the mean value is chosen within la deviation flSheth fe Tormen 
20041 1. for halos with mass M. Finally one can write, 


dM lf avdN^E) M A vir p(z) / , ( , , (c vir r_ 2 ) 


dE 


■( E,M,z) = 


dE Ml 

X. 


dc' vir V{c' mr ) 


[h K ir r- 2 )]‘ 


d‘2 (%min 7 Cvir ^—2) 


(5.14) 

In the above r_ 2 is the ratio between 77 ; and r s where 77 ; is the radius at which the 
effective logarithmic slope —2 that follows from the relation, d/dr (r 2 g(r))\ (- 2 ) = 0. For 


„(- 2 ) _ 

S - 

in Eq. 15.141 can be cast into the form, 


NFW profile, r) = r s . Hence c mr r _ 2 = R V i r /r and the form of integration I n {x m i n , x max ) 


%max) 


dxx g n (x) . 


(5.15) 


Plugging the above equation in Eq. 15.41 the anal ytic form of extragalactic gamma-ray flux 
from DM annihilation can be obtained as flUllio et ah 2002) 


dtp^ 

dE { 


av c 


A f dz a + „, A ) 

„ + ~' h(z) dE 


where the expression for A 2 (z) can be given by, 

,is(z, M)f (is(z, M)) 


A 2 (z) = / dM- 


a(M) 


da 


dM 


A 2 


M 


>,M) 


with 


A m(A, M) — 


A V ir(z) [ 1 , / X l 2 {Xmin,c' vir {z,M)r- 2 ) , , , , \3 

/ ^vir ' VuirJ ,, , , , x -,2 \^vir D A* 7 -* ) 2 ) 

3 J [h(x min ,c' vir {z,M)r- 2 )\ 


(5.16) 


(5.17) 


(5.18) 


In all of the above the concentration parameter, c V i r is defined as 


Rvir 

°vir — ,_ 2 s 

Ts 


(5.19) 


We h ave chosen two forms of co ncentration parameter Cm, r follow i ng M accio et al. flMaccio' et al 
2008) and power law model ( Neto et al. 2007 : Maccio’ et ah 2008). For the first choice 
(Maccio et al), c vir (M,z) = k 200 (R{ z f(M))/'H( y z)) 2 ^ 3 , where k 2 oo — 3.9, 'H(z) = H(z)/H 0 
and z c (M) is the effective redshift during the formation of a halo with mass M. In the power 
law model (second choice) however the expression of c vir (M,z ) is adopted as c vir (M,z) = 
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6.5 'H(z) 2 / 3 (M/M*) 01 , M* = 3.3710 12 h 1 M 0 . This choice of c vir (M,z) provides a rea¬ 
sonable fit within the resolved mass range in the simulations. 


The dark matter substructures are present within halo and form bound objects. The 
mass of the smallest possible such bound object (subhalo) is denoted as M min . The value 
of this minimum subhalo mass, M m ; n is determined from the temperature at which the DM 
particles just start decoupling kinematically from the cosmic background. 


In this w ork we perform our analysis for two typ ical values of M min ; M, 


and 10 9 M 0 ([Martinez et al.ll2009t lBringma/rmll200 
to these subhalos depends inversely on M min . 


r 


= 1O _6 M 0 
The boost factor for y-ray flux due 


5.2. Non-DM Contributions in DIGRB 


The extragalactic gamma-ray spectrum in the energy range between ~ few hundred 
MeV and ~ few hundred GeV as observed by the Fermi-LAT telescope is found to follow 
almost a power law spectrum (^7 oc E~ 2A1 ). There are contributi ons from astrophysical 


sources other than that from possible dark matter annihilation (ITavakoli et al.1120141 1. The 


possible sources that contribute to the diffuse y-ray background other than the DM include 
BL Lacertea objects (BL Lacs), flat spectrum radio quasars (FSRQs), millisecond pulsars 
(MSPs), star forming galaxy (SFG), Fanarof-Riley (FR) radio galaxies of type I (FRI) and 
type II (FRII), ultra high energy cosmic rays (UHECRs), gamma ray bursts (GRBs), star 
burst galaxy (SBG), Ultra High Energy protons in the inter-cluster material (UHEp ICM) 
and gravitationally induced shock waves IGS) etc. The spec tral features of photon spec¬ 
tra originated from the non-DM objects flTavakoli et al.l 1201411 considered in this work are 
concisely summarised in Table [2 


We add up the contributions to EGB both from the annihilating DM in IHDM frame¬ 
work (~ 63.5 GeV DM considered in this work) and the other possible non-DM astrophysical 
sources. The comparison of the sum total value of the y-ray flux with the observed EGB by 
EGRET and Fermi-LAT is shown in Fig. [TO] Needless to mention here that the four plots in 
Fig. [10] are for different parametrisations of concentration parameter c„j r and subhalo mass 
M min . As mentioned earlier we have considered BL Lacs, FSRQs, MSPs, SFGs, FR (type 
I and II), UHECR, GRBs, SBGs, UHEp interacting with ICM and IGS as contributors to 
EGB other than DM and their contributions to EGB are shown with different lines in Fig. [TUI 
The computed total photon spectra in the plots of Fig. [TO] are shown by black solid lines 
while the red solid line is for the minimal non-DM contribution to EGB. The black lines are 
found to be on top of the red lines for Maccio et al. models. 
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We also mention that the extragalactic 7 -ray signal is analysed within the dark matter 
annihilation scenario in Ref. flCholis et al.l 12014 ). In their analysis they have constructed a 
model for the non-DM astrophysical contributions to the extragalactic 7 -ray backgnxmd and 


they have also adopted a 

substructure model based on numerical simulations. The authors 

in Ref. ( 

Cholis et a 

2014 

) have considered subhalo boost factor b s h in their analyses to be 

the following form 

4ndo & Komatsu 

2013) 


b sh (M) « 110 x (M 200 / 10 i2 M 


\0.39 


(5.20) 


where M 2 00 denotes the mass enclosed within a radial region where the avaraged density is 
200 times more than the critical density of the universe. We have performed the calculation 
of extragalactic photon flux for DM for the best-fit model parameter in IHDM based on their 
extragalactic modelling. The result is shown in Fig. [TT] I 11 Fig. [lT]we have only considered 
the contributions to EGB only from radio galaxies, BL lacs, FSRQs and SFGs other than 
that from DM annihilation. The sum total contributions to EGB is found to fit reasonably 
well with Fermi-LAT data. 


5.3. Galactic (sub)halo contribution to DIGRB 


There could be a significant contribution to DIGRB which is of galactic origin along 
the line of sight due to the passing of the signal from extragalactic sources through the 
Milky Way galactic halos and subhalos. From numerical simulation the main DM halo is 


found to host a larg e amount of substructures in form of subhalos flSpringel et al 
Diemand et ah 2007). 


2008b 


The signal from the DM annihilation at the galactic substructures could potentially give 
rise to an almost isotropic signal since this generated 7 -ray flux is propor tional to the less 
centrally concentrated number density distribution of subhalos flSpringel et al.ll2008al ). The 
averaged photon intensity from DM annihilation in such smooth halo of the Milky Way can 
be written as, 


7 f s 1 r 1 ( Ft, 


dE~ 


(aw) dNj 1 
2wip M dE^ VL e 


/ ^ 4^ ( 5 - 21 ) 

i "»DM UIJ 7 “e A, 

where s and D e are the distance from the galactic centre and the observed solid angle. I 11 
the above b and i are galactic coordi nates ( l atitud e and longitude respectively) chosen to 
be 30° < |6| < 90°, 0 < i < 2n f Cholis et ah 2014 ). 77. mw = 21.5 kpc, r viT . .Mw = 258 kpc, 
Ps,mw = 4.9 x 10 6 M 0 kpc -3 , Mvir, mw = 1-0 x 10 12 M q f Klvpin et ah 2002) are chosen for 


































our calculation. 


The photon flux produced in the smooth halo component are much subdominant than 
that yielded in the subhalos and hence they contribute negligibly to extragalactic y-ray back¬ 
ground. In ACDM cosmology, the formation of the structures is assumed to be hierarchical. 
The smaller DM halos are formed first and the larger ones later. In the period of structure 
formation the smaller halos are tidally disrupted after being captured by the larger host 
halos of galaxies and clusters and hence the outer low density layers are stripped in this 
process. Thus the central dense cores only survive and behave as subhalos of the host halos. 
These substructures of DM halo enrich DM phenomenology by giving rise to substantially 
enhancement of the DM annihilation rates within a halo. The contribution to differential 
gamma-ray flux from subhalo can be obtained from the differential luminosity profile of each 
subhalo which is given by, 


dLry 

dE~, 


(aw) dN 1 

2 m DM dEry 


dVshp 2 sh - 


(5.22) 


For an individual subhalo with mass M, the photon intensity can be written as, 


dZ(E 7 , s, M) 1 dL(Ey, M ) 
dE 7 47 ts 2 dE y 


1 b gs (av) dNy M 2 
4ns 2 2 m^ M dEy r SjSh (M) 3 


S[Ccut(M)]. 


(5.23) 


where r S;S h denotes the scale radius of the subhalo 
the contribution from substructure within each subhalo 


In the above the factor b gs determines 
(‘subsubhalo’) and is chosen to be 


2 ( Kuhlen et ali 2008). The function Q[c cut (M)\ which can be obtained using integral over 
the volume of each satellite and the form of subhalo concentration c cut 0 following Ref. ( Ando 


20091 ). can be given as, 


G[Ccnt(M)} 


1 

127T 


1 

(1 + C cu t) . 


ln(l T Ccut) 


Ccut 

1 T Ccut 


-2 


(5.25) 


The total y-ray intensity at Earth from the annihilation of dark matter particles in 
galactic subhalos can be written after integrating Eq. 15.231 over the distribution of Milky 


^The NFW halo profile of Milky Way is chosen as p M wW = ( r /r s M w)(m/r. M w) where r and r s,MW are 
the galactocentric and scale radii respectively and p Sj mw is the scale density. The above-mentioned halo 
profile is assumed to extend up to the virial radius r v ; rj MW with the virial mass M v i r ,MW- 

** We choose the DM density within each subhalo of mass M to be truncated NFW halo profile, 

/ IA / PNFw(r s h|M) for r sh < r cut , . . 

PshRsh M) = < . (5.24) 

l 0 for r sh > r cu t, 


where c cut is the cutoff radius for this profile. 























21 


Way subhalos as, 

c?X sh (E 7 ) _ r dn&{M, s, £, b ) dl(E 1 , s, M ) 

dEry J dM dEry 


(5.26) 


where J dMdV(dn s h/dM) is the total number of subhalos in the Milky Way. The form of 
subh alo mass function, dn^/dM is chosen to be the anti-biased model 0 following Ref. ( Ando 


20091 ) for our calculation. 


In order to confront observations, we are interested in the averaged intensity of y-rays 
per unit energy emitted due to DM annihilation over the whole galaxy, 


dd K h ( Er-y ) 
dE~y 


1 dl sh (E 7 

dE~ 


1 

a 


dV dM 


dn s h(M, s, b ) dZ(E 7 , s, M) 


e J M* JV*(M) 


dM 


dV 


dn s h(M, s,£,b) 1 


dM 

(crv) dN 1 


dE~ 


M 2 


I M* 




dM Ati s 2 m() M r S;S h(M) : 


r ^[c cut (M)|5.27) 


where K is the volume beyond which satellites remain unresolved. The considered mass 
range of the subhalos is 1O _6 M 0 < M* < 1O 1O M 0 . Since from the luminosity L one gets the 
knowledge of the subhalo mass M we consider the subhalo mass range in such a way that the 
bright as well as the faint subhalos are included in the calculation. Also since luminosity is 
directly related to the flux sensitivity of Fermi (F sens ) by the relation, L(M) = 47rs^(M)F sens , 
they remain unresolvable beyond the distance s*(M) = \/'L(M) /47rF sens , where the flux 
sensitivity of Fermi-LAT, F sens = 2 x 10 -lo cm -2 s -1 (Ando! 2 0091 ) and L(M ) is the luminosity 
obtained by integrating Eq. 15.221 over energy. 


In Fig. [T21 we show the contributions to the EGB from both galactic smooth halos 
and subhalos. From the left panel of Fig. [12] we see that the contribution to the EGB is 
much subdominant for the annihilation of DM within galactic smooth halo. For our chosen 
galactic substructure model, this contribution is comparable to that from extragalactic dark 
matter annihilations calculated using c V i r of Maccio et al. model. For the other case where 
the gamma ray flux from extragalactic dark matter annihilations have been calculated using 
c V i r of power law model, the contributions from the galactic subhalos are much negligible. 


'Tn another model (‘unbiased’) for n s h{r), the subhalo distribution is assumed to follow its parent NFW 
halo distribution whereas in the anti-biased model the subhalo distribution is flatter than NFW halo (Ando 


2009). 
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6. Summary and Conclusion 

We have chosen a simple dark matter (DM) model, namely inert Higgs doublet model 
(IHDM) where scalar sector of Standard Model is extended by adding another SU(2 )l dou¬ 
blet. The newly added doublet does not generate any VEV after spontaneous symmetry 
breaking. The ‘inert’ doublet is considered to be the DM candidate. The stability of DM is 
ensured by imposing discrete Z 2 symmetry. The model, in general, provides a broad range 
of DM mass from GeV to TeV range. In this study we only consider the lower mass range of 
DM in this model. The analysis of experimental data for DM relic density from PLANCK 
experiment and the other direct detection experimental results for the case of this IHDM 
gives a set of best fit values for DM mass, annihilation cross section and other model pa¬ 
rameters. We adopt this best fit point (obtained using y 2 minimisation) for IHDM mass 
from such analyses. Thus the DM mass of ~ 63.5 GeV is our chosen benchmark point in the 
present work. We study the 7 -ray spectrum obtained from the annihilation of this chosen 
DM particle in IHDM framework and interpret various types of continuum 7 -ray fluxes with 
astrophysical origins measured by Fermi-LAT satellite. 

In this study we compare our calculated 7 -ray flux with the galactic centre 7 -ray excess 
in the light of this model. For this we have employed different analysed Fermi-LAT residual 
7 -ray flux data for different angular regions around the galactic centre. The calculated 
low energetic photon spectra from the annihilation of the DM particle with benchmark 
value of mass (~ 63.5 GeV) in IHDM for various chosen regions surrounding the galactic 
centre are found to be in the same ballpark as reported by these studies. Although in some 
previous analyses it was argued that the photon spectra originated from different annihilation 
channels of dark matter particles with low masses can possibly fit the obtained data, very 
recent analyses have obtained the resulting best fit masses of dark matter to be much more 
conservative (and also somewhat higher as well). We have computed the photon spectra for 
our benchmark scenario in IHDM framework and have confronted with the residual photon 
spectra obtained for all of the above-mentioned studies. Our theoretical calculations for 
photon spectra in this model have been performed after suitable parametrisation of the dark 
matter halo parameters, region of interest surrounding the galactic centre etc. 

We also address the prospects of the continuum 7 -ray signal which may come from 
DM-dominated dwarf spheroidal galaxies (dSphs) in case they originate from dark matter 
annihilation. We then compare the 7 -ray flux that can be obtained from the benchmark 
IHDM dark matter with mass ~ 63.5 GeV. For this we choose 18 Milky Way dSphs whose 
J-factor can be estimated from measurements. The uncertainties in the measurement of 
J-factor for different dSphs are also incorporated in our calculations. The calculated photon 
spectra for IHDM benchmark point are seen to obey the allowed limits for observed spectra 
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of continuum 7 -ray. 

After addressing the issues regarding indirect DM searches with 7 -ray signals from 
various galactic cases, we finally confront the extragalactic 7 -ray signal with that from the 
annihilation of low mass DM (considered in this work) in IHDM scenario. We calculate 
the extragalactic 7 -ray flux for different extragalactic parametrisations and compare with 
the observed extragalactic gamma ray background by EGRET and Fermi-LAT. For this 
we consider several possible classes of non-DM astrophysical sources which may yield 7 -ray 
signal embedded in the extragalactic background. Although there are too many uncertainties 
involved in modelling of such astrophysical sources and other parameters for extragalactic 
flux calculation, we have shown that the considered low mass DM in IHDM can generate 
photon flux within the observed flux limit. 

From the detailed study of various galactic and extragalactic 7 -ray searches for probing 
the indirect signatures of DM in light of IHDM, we can conclude that the low mass DM 
in this model framework is still a viable candidate to be probed in future 7 -ray searches. 
Although we have performed the thorough analysis considering only a single Higgs-portal 
model, the analysis is valid for any simple Higgs-portal DM model such as singlet scalar 
DM model, singlet fermion DM, inert Higgs triplet model etc. with DM mass in the same 
ballpark as in our study. 
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Table 1: Benchmark Point : IHDM Parameters 


Model 

M h o (GeV) 

M h ± (GeV) 

M h o (GeV) 

M a o (GeV) 

Al 

A 2 

Parameters 

126.016 

73.78 

63.54 

166.16 

-3.29 x 10~ 3 

5.67 x 10~ 4 


Table 2: Benchmark Point : Observables 


DM 

Q.h 2 

(av) (cm 3 s *) 

csi 

(Pb) 

Observables 

0.1173 

2.37 x 10 -26 

8.89 x 10 _u 

Annihilation 

H°H° bb 

H°H° W+W~ 

H°H° gg 

H°H° ll 

H°H° ->• cc 

H°H° ZZ 

Cross-section 

69.2% 

9.61% 

9.49% 

7.37% 

3.29% 

0.48% 


This preprint was prepared with the AAS IATpX macros v5.2. 
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Table 3: Limits on DM annihilation cross-section from 7 -ray flux limits for various dwarf 
spheroidal galaxies for the benchmark DM mass of Table [T| 


dSphs name 

longitude 
l (deg) 

latitude 
b (deg) 

distance 

(kpc) 

D ffl0 (J NFW ) 

(Zogi 0 [GeV 2 cm -5 sr]) 

log, o(aNFW) 
(D(/i 0 [deg]) 

upper limit on 
(ov) (cm 3 s -1 ) 

Bootes I 

(DalhOra et ah 20061 

358.1 

69.6 

66 

18.8 ±0.22 

-0.6 ±0.3 

2.33 x 10“ 24 

Bootes II 

353.7 

68.9 

42 

- 

- 

- 

Bootes III 

35.4 

75.4 

47 

- 

- 

- 

Canes Venatici I 
('Simon & Geha 2007) 

74.3 

79.8 

218 

17.7 ±0.26 

-1.3 ±0.2 

9.65 x 10 -25 

Canes Venatici II 
('Simon & Geha 2007) 

113.6 

82.7 

160 

17.9 ±0.25 

-1.1 ±0.4 

8.14 x 10 -25 

Canis Major 
('Walker et ah 2000a,) 

240.0 

-8.0 

7 




Carina 

('Simon & Geha 2007) 

260.1 

-22.2 

105 

18.1 ±0.23 

-1.0 ±0.3 

2.28 x 10“ 25 

Coma Berenices 
('Simon V Geha 2007) 

241.9 

83.6 

44 

19.0 ±0.25 

-0.6 ±0.5 

1.11 x 10“ 24 

Draco 

(Munoz et ah 2005) 

86.4 

34.7 

76 

18.8 ±0.16 

-0.6 ±0.2 

3.87 x 10 -25 

Fornax 

(Walker et ah 2000a,) 

237.1 

-65.7 

147 

18.2 ±0.21 

-0.8 ±0.2 

2.53 x 10 -25 

Hercules 

('Simon & Geha 2007) 

28.7 

36.9 

132 

18.1 ±0.25 

-1.1 ±0.4 

9.97 x 10“ 26 

Leo I 

(Mateo et al. 2008) 

226.0 

49.1 

254 

17.7 ±0.18 

-1.1 ±0.3 

4.37 x 10“ 25 

Leo II 

(Koch et ah 2007) 

220.2 

67.2 

233 

17.6 ±0.18 

-1.1 ±0.5 

3.88 x 10 -25 

Leo IV 

('Simon & Geha 2007) 

265.4 

56.5 

154 

17.9 ±0.28 

-1.1 ±0.4 

3.72 x 10 -24 

Leo V 

261.9 

58.5 

178 

- 

- 

- 

Pisces II 

79.2 

-47.1 

182 


- 

- 

Sagittarius 

5.6 

-14.2 

26 

- 

- 

- 

Sculptor 

('Walker et ah 2000a,) 

287.5 

-83.2 

86 

18.6 ±0.18 

-0.6 ±0.3 

3.41 x 10“ 24 

Segue 1 

('Simon et ah 2011) 

220.5 

50.4 

23 

19.5 ±0.29 

-0.4 ±0.5 

1.16 x 10“ 24 

Segue 2 

149.4 

-38.1 
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Sextans 

(Walker et, ah 2000a,) 

243.5 

42.3 

86 

18.4 ±0.27 

-0.9 ±0.2 

1.14 x 10 -25 

Ursa Major I 
('Simon & Geha 2007) 

159.4 

54.4 

97 

18.3 ±0.24 

-1.0 ±0.3 

1.64 x 10 -25 

Ursa Major II 
('Simon & Geha 2007) 

152.5 

37.4 

32 

19.3 ±0.28 

-0.5 ±0.4 

1.33 x 10 -24 

Ursa Minor 

105.0 

44.8 

76 

18.8 ±0.19 

-0.5 ±0.2 

6.54 x 10 -24 

Willman 1 

(Willman et ah 2011) 

158.6 

56.8 

38 

19.1 ±0.31 

-0.6 ±0.5 

4.03 x 10 -24 
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Table 4: Overv iew o f the minimal non-DM contributions to the total extragalactic 7 -ray 
background (Ta vak oli et ah 2014 ) 


Non-DM objects 

Photon Energy Spectra in GeV 1 cm 2 s for ] ) 

BL Lacs 

3.9 x 10~ 8 E- 2 ' 23 

FSRQ 

3.1 x 10- 8 ^- 2 ' 45 

MSP 

1.8 x 10- 7 E- 1 - 5 exp (-§) 

SFG 

1.3 x 10" 7 E- 2 ' 75 

FR I & FR II 

5.7 x 10- 8 E- 2 - 39 exp(-|^) 

UHECR 

4.8xlO- 9 E- 4 - 8 exp - (Jfo) 

GRB 

8.9 x lO- 9 ^- 2 - 1 

SBG 

0.3 x 10~ 7 E~ 2A 

UHEp ICM 

3.1 x 10~ 9 £- 2 ' 75 

IGS 

I fel 2 04 for E L < 10 GeV 

0.87 x 10~ 10 x < 2 . 13 \ 

(■for E 7 > 10 GeV 
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Fig. 1.— Various DM direct detection experimental bounds on the benchmark point of Ta- 
ble Ul See text for details. 



Fig. 2.— Residual GeV gamma-ray flux from the inner 5° surrounding the galactic centre. 
See text for details. 











35 



Fig. 3.— Left panel: Comparison of calculated 7 -ray flux (green line; assuming DM annihi¬ 
lation) with the observed ~ 1 — 3 GeV gamma-ray excess from GC (red data points). Right 
panel: Same as the left panel but here the observed data are for galactic ‘North-South’ region 
(\b | < \(\; red data points) and for galactic ‘East-West’ region (\b\ > \£\; blue data points). 
See text for details. 




Fig. 4.— Left panel: Comparison of calculated 7 -ray flux (red points) with the observed 
residual 7 -ray spectrum for the region in galactic coordinate \i\ < 20°, 2° < \b\ < 20°. The 
black line is for canonical J-factor while green and blue lines are calculated with minimum 
and maximum deviations of J-factor from its canonical value. Right panel: Comparison of 
the calculated results in the left panel with observed 7 -ray spectra (red points) obtained by 
studying 15° x 15° region around GC. See text for details. 
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Fig. 5.— Comparison of computed j-ray flux with the bin-by-bin integrated 7 - ray energy-flux upper limits for each 
dSph. The downward red-coloured arrows represent bin-by-bin upper limits on the 7 - ray energy-flux at 95% CL. The 
blue lines denote the j-ray fluxes calculated using the central values of integrated J-factor whereas the green band is for 
the uncertainties in the measurement of integrated J-factors for 18 Milky Way dSphs. See text for details. 
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Fig. 6.— The variation of optical depth (transparency coefficient) e T with the energy at 
detection (E) and the redshift (z) of photon emission for minimum UV background model. 
The black zone in this figure represents total opaque region while the yellow zone is for total 
transparent one. See text for details. 



Mass (M sun hi 1 ) 



Fig. 7.— The variation of the linear power spectrum P(k ) of matter density perturba¬ 
tions with the wavenumber k of the fluctuations for different redshifts is shown in the 
left panel. In the right panel the variance a of the density perturbations is shown as 
a function of halo mass for different redshifts. In both plots the values of redshift z = 
10" 3 ,10~ 2 ,10 -1 ,10°, 10 1 ,10 2 and 10 3 . See text for details. 
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Fig. 8.— The fraction of mass collapsed, f(a) in Sheth-Torman model for different redshifts 
z and the halo masses M is shown in left panel. The variation of Sheth-Torman halo mass 
function ^ with the redshift z and the halo mass M is shown in right panel. See text for 
details. 



Fig. 9.— Comparison of the computed extragalactic gamma-ray fluxes from the annihilation 
of ~ 63.5 GeV dark matter (benchmark point) in IHDM framework for different extragalac¬ 
tic parametrisations. We consider two models of the concentration parameter c V i r , namely 
a)power law model and b)Maccid et al. model. Also the minimum extragalactic subhalo 
mass M min are chosen to be 1O~ 6 M 0 and 10~ 9 M o . Calculation with power law model yields 
enhanced y-ray flux compared to that with Maccio et al. model. For low M min , y-ray flux 
increases and this enhancement is smaller for Maccio et al. model compared to that for power 
law model. See text for details. 
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c vh . = Power Law Model, M ■ = 10 M 


c • = Power Law Model, M ■ =10' M 
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Extragalactic y-ray Background measured by Fermi-LAT - 

Extragalactic y-ray Background measured by EGRET - 

BL Lacertae Objects - 

Flat Spectrum Radio Quasars 
Millisecond Pulsars 
Star Forming Galaxies 

Fanarof-Riley I & II Radio Galaxies . 

Ultra High Energy Cosmic Rays . 

Gamma Ray Bursts - 

Star Burst Galaxies . 

Interaction of Ultra High Energy Protons with Inter-Cluster Material ---- 

Gravitationally Induced Shock Waves in the Intergalactic Medium - 

Total Minimal Non-DM Contribution to Extragalactic Background - 

y-ray from Annihilation of63.54 GeV DM in IHDM Framework . 

Total Contribution (DM + Non-DM) to Extragalactic y-ray Background - 


Fig. 10.— The observed extragalactic y-ray fluxes by EGRET and Fermi-LAT are compared 
with the sum total value of the y-ray fluxes obtained from the present calculation. The calcu¬ 
lated value of extragalactic y-ray flux is obtained by summing over the y-ray flux calculated 
from DM annihilation for IHDM LIP DM (considered in this work) and other possible y- 
rays (extragalactic y-ray sources of non-DM origin) from extragalactic sources. See text for 
details. 
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Fig. 11.— Comparison of the observed 7 - ray fluxes by EGRET and Fermi LAT theoret¬ 
ical results for 7 - ray spectra for ~ 63.5 GeV DM in IHDM consider i ng the modelling of 
extragalactic and astrophysical parameters as done in Ref. Gholis et al. (201$) . See text for 
details 



E y (GeV) 



Ey(GeV) 


Fig. 12.— Galactic smooth halo and subhalo contributions to the extragalactic gamma-ray 
flux for LIP ~ 63.5 GeV dark matter in IHDM. See text for details. 



























